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Photon-based quantum information processing is based on manipulating multi photon interference. We 
focus on the Hong-Ou-Mandel (HOM) dip in the photon coincidence rate which provides a direct measure 
of interference of indistinguishable photons linked to their Bose statistics. The effect has been first observed 
with entangled photons generated by parametric down conversion and then extended to independent 
emitters. Fluctuations caused by coupling between emitters and a bath can erode the interference which 
causes the dip. Here we show how the magnitude and shape of the HOM dip is affected by spectral diffusion 
induced by coupling to a brownian oscillator bath. Conditions for maintaining and controlling the 
interference are specified. 

The ultrafast optical control of quantum interference of single photons is of fundamental interest with broad 
applications to entanglement generation 1 , quantum computing 2 " 4 , communication 5 , and information pro- 
cessing 6-8 . The HOM dip in the photon coincidence counting (PCC) rate 1 is an ultimate measure of photon 
indistinguishability 910 . It has been observed in a broad variety of systems including remote frequency- detuned 
organic molecules 11,12 , quantum wells 13 , quantum dots 14 " 16 , nano crystals 17 , nitrogen vacancies in diamond 18 " 20 
and carbon nanotubes 21 . We examine photon interference in the setup shown in Fig. la. A pair of photons is 
generated by two remote two-level molecules a and b with ground^ and excited state e a , oc = a,b. These photons 
then enter a beam splitter and are subsequently registered by time-and-frequency gated detectors s and r. There 
are two types of possible outcomes: two photons registered in one detector (s or r) or coincidence where one 
photon is detected in each. The ratio between these outcomes reflects the Bose statistics and degree of distinguish - 
ability of photons. If the two photons incident on the beam splitter are indistinguishable the PCC vanishes. This is 
known as the HOM dip. The dip is displayed by varying the position of the beam splitter which causes delay T 
between the two photons. The normalized PCC rate varies between 1 for completely distinguishable photons and 
0 when they are totally indistinguishable. For classical fields and 50 : 50 beam splitter the PCC rate may not be less 
than 1/2. We denote the photons to be indistinguishable (distinguishable) if the PCC rate is smaller (larger) than 
1/2. 

PCC is typically measured using pure time-domain detection 22 ' 23 . Originally performed with entangled 
photons generated by parametric down conversion (PDC) 1 the shape of the dip vs delay is usually related to 
the two-photon state envelope which is governed by an effective PDC Hamiltonian 24 . Bath induced fluctuations 
can become important for remote emitters and have been introduced phenomenologically 12 . Here we present a 
microscopic theory for PCC with bath fluctuations by formulating the signal in the joint field-matter space. 

Results 

Generation and detection of PCC by two remote emitters. The time-and-frequency gated PCC signal is 
described by the two pairs of loop diagrams shown in Fig. lc. Each loop represents molecule (a or b) which 
undergoes four field-matter interactions and each detector interacts twice with the field. Fig. lc shows that after 
interacting with pump (with its ket) at time t 2 molecule a evolves in the coherence p 6aga during time interval t' 2 . 
The second interaction of the pump with the bra then brings the molecule into a population state p^ which 
evolves during interval t Y until the first interaction with spontaneous emission mode occurs with ket. The 
molecule then evolves into the coherence p during t \ until the second bra- interaction of spontaneous 
mode. During population and coherence periods, the characteristic timescale of the dynamics is governed by 
population relaxation and dephasing, respectively. 
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Figure 1 | Time-and-frequency resolved measurement of PCC with spectral diffusion. Schematic of the PCC experiment with two source molecules - 
(a), the two-level model of the molecule with SD used in our simulations - (b). (c) - Loop diagrams for the PCC rate of emitted photons from two 
molecules (for diagram rules see 26 ). The left and right branches of each diagram represent interactions with ket- and bra- of the density matrix, 
respectively. Field-matter interactions with the pump pulses p a and pb (blue), spontaneously emitted s, s', r, r' photons (red) and detectors (brown). 



All relevant single-molecule information is given by the four point dipole 
correlation function F a (t u t 2 ,t 3 ,t 4 ,) = <V^(fi) V^(f 2 )V^(f 3 ) V^(f 4 )) a , 
where V and V 4 are the lowering and raising dipole transition opera- 
tors, respectively. Diagrams i in Fig. lc represent non-interfering 
term given by a product of two independent fluorescence contribu- 
tions of the individual molecules. Diagrams ii represent interference 
described in the joint space of the two molecules and involve the 
interference of eight quantum pathways (four with the bra and four 
with the ket) with different time orderings. Each molecule creates a 
coherence in the field between states with zero and one photon |0)(1 1 
and 1 1)(0|. By combining the contributions from a pair of molecules 
we obtain a photon population |1)(1| that can be detected 24 ' 25 . For a 
pair of identical molecules, the beam splitter destroys the pathway 
information making the molecules indistinguishable and giving rise 
to quantum interference. 

Bath induced fluctuations via e.g. spectral diffusion (SD) can des- 
troy the quantum interference. Generally SD can manifest itself as 
either discrete random jumps of the emission frequency 21 ' 31,32 or as a 



broadening of a hole burnt in the spectrum by a narrowband 
pulse 33 ' 34 . We focus on the SD in the "hole burning" limit (HBL) 
(see Methods for the precise conditions and Section S3 of SI). This 
is relevant to the crystals which store information in the form of 
reversible notches that are created in their optical absorption spectra 
at specific frequencies. Long storage times 35 , high efficiencies 36 , and 
many photon qubits in each crystal 37 can be achieved in this limit. 

We have calculated the PCC signal (see Methods) for the output 
fields £ 3 and E 4 of the beam splitter (see Fig. la) that contains all 
relevant field matter interactions (see Section S3 of SI) using the SD in 
HBL. It first assumes that the dephasing is much faster than the 
fluctuation timescale, i.e. f^^cA" 1 , k= 1, 2, 3, 4. Second, if excitation 
pulse duration o~ l and the inverse spectral (crjj \ and temporal 

(°r) >j = r? s £ ate bandwidths of the detectors and much shorter 

than the fluctuation time scales, one may neglect the dynamics during 
the delay between population evolution and its detection. Under these 
conditions the PCC signal is given by 
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Figure 2 | Effect of gating spectrograms on photon indistinguishability in the absence of spectral diffusion. PCC signal for A a = A b = 0 given by 
Eq. (4) vs the delay between the detection t s = t r for different values of the time - (a) and frequency - (b) gating bandwidths for the fixed frequency 



crjy = 1 MHz - (a) and time 0*^ = 1 MHz - (b) bandwidths, j = r, s. The parameters are chosen as = 3.2 GHz, co\ 



: 10 MHz, T= 10 ns, o p 



3 MHz, A = A h = 25 MHz, co v 



-Ah 
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(1) 



where expressions in the last line represent permutation of the mole- 
cules a and b, Tj = {tj,cdj} represents a set of gating parameters for 
the detector j = r, s. C a (r = {f, w}) is the time-and-frequency 
resolved fluorescence of molecule oc = a,b corresponding to diagram 
i in Fig. lc: 



d t (t,<o) = d a0 (t)e 



-la) (gj-co^t)) 2 



(2) 



CL^ = <x> a — 2 a is the mean absorption and fluorescence frequen- 
cy. l{(Ti,t 2 ,T) and l{(ti,T 2 ^) with t x < t 2 are the interference 
contribution a = a, b, j = r, s corresponding to diagram n in 
Fig. lc 



4(r 1? f 2? T)=4(M 2 )e 



4< 



|-K«(^)V 



x e 



(3) 



l/(r r ,r s ;T)=© B (? 5 -f r )+c^ fl (f r J I ,© r )T+(^/A a )(2^ 

A 2 

+ Fa&-tr))-{a++b 9 r++s), ^( t )=J2 0C=a ^ F M with 



^a(0 =e ~ a +A a f— 1, a = a, ^ and all the remaining parameters 
are listed in Eqs. (S24)-(S31) of the SI. The contribution of Eq. (2) 
enters signal in Eq. (1) as the amplitude square coming from each 
molecule in the presence of fluctuations. The interference term (3) 
generally cannot be recast as a product of two amplitudes 26 . In the 
following we simulate the results of Eqs. (l)-(3) using the typical 
parameters of the TPI experiments 9 " 20 . 

Variation of the HOM dip with gating. We first examine the effect 
of time and frequency gating on photon indistinguishability and 
PCC. In the absence of the SD and using identical detectors 
(T r w = a s w , <7 r T = a s T the signal (1) for tj = t s — t r is reduced to 



£ c 34 (t; T) - 1 - y\ cos(w ab t d + Q T T)e~ 



-K A «+ A2 0^ 5 (4) 



where the parameters r], Q T and o x are given in Eqs. (S32)-(S35) of 
the SI. Eq. (4) is commonly derived in the field space alone by adding 
a phenomenological Gaussian function of the time delay 1 . In the 
present analysis the Gaussian shape is characteristic to the 
harmonic bath. 

Fig. 2a depicts Eq. (4) vs the delay between detectors t s — t r for 
different values of the time gate bandwidth a T and fixed spectral gate 
bandwidth <j oj . If the time gate bandwidth is smaller than the splitting 
w° b — w° a = 10 MHz the HOM dip is shallow and does not go below 1/ 
2 which implies that photons generated by the two molecules are 
distinguishable. If the time gate bandwidth is larger than the splitting, 
the HOM dip drops below 1/2 indicating that the photons are indis- 
tinguishable. Similarly in Fig. 2b we fix the time gate bandwidth a T 
and vary the frequency gate o OJ . As o OJ is increased the photons 
gradually become more indistinguishable and the dip becomes lar- 
ger. In the original HOM experiment 1 the PCC was plotted vs the 
delay T. Here we fix T and vary the delay between the detectors 
td = t s — t r . In the absence of fluctuations both T and enter the 
signal (4) in a similar way. However once fluctuations are included 
Eq. (1), t d is directly connected to fluctuations whereas the delay Tis 
not. Below we study the PCC variation with t d holding T fixed. 

Narrowband spectral filters are typically used in the experiments 
with broadband laser pulses. The broadband excitations degrade the 
quantum interference 27 " 29 . This can be easily explained since the in 
addition to selecting a well defined frequency, the narrow frequency 
gate makes the ultrashort pulse longer and therefore two such 
photons are more likely to overlap in time and harder to distinguish, 
enhancing the HOM dip. In the present application we use narrow- 
band excitation pulse and the degree of distinguishability is con- 
trolled by the interplay of spectral and temporal gating widths 
through the Wigner function. Therefore if the time gate is narrow, 
the broad frequency gate makes photons less distinguishable increas- 
ing the visibility of HOM dip. 

Time-and-frequency resolved fluorescence with spectral dif- 
fusion. The simplest way to observe SD is by time-and-frequency 
resolved fluorescence. The molecular transition frequency is coupled 
linearly to an overdamped Brownian oscillator that represents the 
bath (see Fig. lb). The absorption and emission lineshape functions 
for a pair of molecules are given by 30 
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Figure 3 | Effect of SD on individual molecule emission lineshapes. 

Absorption (5) - (a) and fluorescence (6) - (b) line shapes vs displaced 

frequency Acq = co — - (co a + co^). Time-and-frequency resolved 

fluorescence (2) displayed as a snapshot spectra for molecule a - (c) and 
b - (d) and depicted as a two-dimensional contour plot - (e) and (f), 
respectively. Parameters of the system are: co° a = 3.2 GHz, co 0 b — co 0 a = 
600 MHz, A fl = 220 MHz, A b = 205 MHz, A a = 5 MHz, A b = 3 MHz, 
co p = w° b + 4, o-p = 102 MHz, d T = a j a = 100 MHz, j = r, s, and the 
temperature is set to 10 K. 



a A {co)= (2nA,r 1/2 e~ 



= a,b 



<x = a,b 



(6) 



where 2k a is the Stokes shift and A a = \/2X a k B T /H is a linewidth 
parameter. Together with the relaxation rate A a (see Eq. (11)) these 
parameters completely describe the SD model and govern the 
evolution of the emission linewidth between the initial time given 
by Eq. (5) (see Fig. 3a) and long time given by Eq. (6) (see Fig. 3b). 
Eqs. (5) -(6) are obtained in the slow nuclear dynamics limit: 
A a <A a . 

The time-and-frequency resolved fluorescence (Eq. (2)) is 
depicted as a series of the snapshot spectra at different times for 
molecule a in Fig. 3c. It shows a time dependent frequency redshift 
cb a (t) and time dependent spectral broadening given by a a0 (t) (see 
Section S4 of SI). Initially w a (0) =w° a -\- l a whereas at long times 
&>a(°o) = co a ~ K-> where 2k a is the Stokes shift. 2D representation 



of Eq. (2) is given in Fig. 3e. Fig. 3d and 3f show the same signals for 
molecule b. Because of the different reorganization energies X a , k b 
and relaxation rates A a , A b the Stokes shift dynamics and dispersion 
are different. Even when the absorption frequencies are the same co a 
= co b , the fluorescence can show a different profile due to SD. This 
affects the distinguishability of the emitted photons as will be demon- 
strated below. 

Combined effect of gating and SD on the HOM dip. We depict the 
2D fluorescence (2) in Fig. 4a. The vertical line marks the pump 
frequency co p = co° b - tuned midway between absorption and 
fluorescence at long time (see Eqs. (5-6)), the Stokes shift is 
~ 1 kHz. Fig. 4b shows the same for co p = co° b + h where we have a 
larger Stokes shift ~2 kHz. At co p = a)° b + 2k b the Stokes shift 
~3 kHz as shown in Fig. 4c. Fig. 4d-f show the PCC signal (1) for 
the parameters corresponding to Fig. 4a-c, respectively. The HOM 
dip is less pronounced as the degree of indistinguishability drops 
from 0.275 - Fig. 4d, to 0.3 - Fig. 4e further to 0.375 in Fig. 4f. 

The distinguishability is also affected by the molecular transition 
frequencies. Fig. 4g shows that for fixed time and frequency gate 
bandwidths er^ = cr^ = 100 MHz, j = r, s the photons are distingui- 
shabe as long as co° b — co° a > 1 10 MHz and are indistinguishable other- 
wise. The effect of the time delay T caused by the position of the beam 
splitter is depicted in Fig. 4h. For T < 10 ns (> 100 MHz bandwidth) 
the photons are indistinguishable, after that they become distin- 
guishable. The SD timescale is one the key parameters affecting the 
degree of indistinguishability. Using Eq. (12) we fixed the absorption 
linewidth T a and varied A a and A a . The PCC signal (1) depicted in 
Fig. 4i shows that if the molecules have nearly degenerate transition 
frequencies co° b — co° a = 1 MHz for slower fluctuations A a = 0.5 MHz, 
A b = 1 MHz the photons are indistinguishable. Increasing the SD 
rate of one of the molecules increases the degree of distinguishability, 
e.g. for A a = 10 MHz the HOM dip becomes 0.225, which is still less 
than classical 1/2 limit. However for fast SD A a = 18 MHz, the HOM 
dip becomes 0.625 and photons are clearly distinguishable even 
though both time and frequency gates are broader than the difference 
in transition frequencies. 

We further illustrate the effect of frequency and time gating in the 
presence of spectral diffusion. Fig. 4j shows that if two molecules have 
different SD timescales (A a =15 MHz and A b = 1 MHz) and the 
frequency gate bandwidth is narrow (cr^ = 50 MHz) the photons are 
rendered distinguishable and HOM dip is 0.6. By increasing the the 
photons become indistinguishable and HOM dip is 0.48 for o OJ = 
120 MHz and 0.35 for o w = 200 MHz. In all three cases we kept the 
time gate fixed at o T = 100 MHz. Alternatively we change the time gate 
bandwidth while keeping the frequency gate fixed at o OJ = 100 MHz. 
Fig. 4k shows that initially indistinguishable photons at o T — 80 MHz 
with HOM dip 0.675 become indistinguishable at o T = 110 MHz with 
HOM dip 0.45 and at o T = 150 MHz with HOM dip 0.275. Thus, if the 
presence of the bath erodes the HOM dip the manipulation of the 
detection gating allows to preserve the quantum interference. 

Discussion 

To put our ideas into more practical perspective and connect with 
quantum processing technologies, we note that recent progress in 
long-term quantum memories using HBL where entanglement is 
achieved with telecom photons, proved the possibility of quantum 
internet 38,39 . Simultaneous time and frequency detection can be 
achieved using high-speed photodiode which converts fast optical 
signal into a fast electrical signal, fast oscilloscopes to observe the 
waveform, wide bandwidth spectrum analyzers and other elements. 
Short pulse characterization using time-frequency map such as fre- 
quency-resolved optical gating (FROG) 40 , spectral phase interfero- 
metry for direct-field reconstruction (SPIDER) 41 are well established 
tools for ultrafast metrology 42 ' 43 . Extending these techniques to a 
single photon time and frequency resolved detection is challenging 
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Figure 4 | The combined effect of SD and gating on photon indistinguishability. Time-and-frequency resolved fluorescence from two molecules 
(vertical black line represents excitation frequency) - top row for different excitation frequencies: co p = co b - (a), aP h + k h - (b), (o\ + 2k h - (c). The 
corresponding PCC signal vs delay between detector clicks: l s — t r - (d) - (f) (middle row) for the values of the excitation frequencies in (a) - (c), 
respectively. Here co° b — co° a = 100 MHz, T = 1 ns, A a = 5 MHz, A b = 3 MHz. PCC for different transition energies of molecules excited at co p = co° b + k\, - 
(g); PCC for different values of the delay Twith aP b — w° a = 100 MHz - (h); PCC for different values of the SD time scale A a and and fixed linewidth T ay 
T b according to Eq. ( 12) - (i). Excitation pulse bandwidth o p = 102 MHz, gating bandwidths = = 100 MHz, j = r, s, transition energies coj — aP a = 1. 
PCC for different frequency gate bandwidths - (j) at fixed time gate bandwidth a T = 100 MHz and for different time gate bandwidths - (k) at fixed 
frequency gate bandwidth o OJ = 100 MHz. Molecules have distinct SD timescales A a = 15 MHz, A b = 1 MHz and co° b — co° a = 1 MHz. 



and can be achieved if combined with on-chip tunable detectors 44 or 
upconversion processes 45,46 . 

In summary we have employed a time- ordered superoperator 
Green's function formalism to calculate the PCC from two remote 
emitters coupled to a bath. In the absence of a bath, time gating is 
sufficient to reveal quantum interference. Even if the emitters have 
well resolved transition frequencies with splitting Aco quantum inter- 
ference manifests if the temporal resolution At is higher than ItiIAoj. 
Once the bath is included, time-and-frequency gating is necessary to 
resolve the dip. We showed how various model parameters affect the 
indistinguishability. 

Methods 

Coincidence measurements. The PCC signal is defined via the four point correlation 
function of the time-and-frequency gated electromagnetic fields E 3 and E 4 incident on 
the detector 22 

^ 4 (r r ,r s ) = 

r dt s r ^(^w^w^w^w), (7) 

J — oo J — oo ^ ' 

where (...) =tr(. .. p) is tracing with respect to p - the total field plus matter density 
operator. The fields in the output 3, 4 and input 1, 2 ports of the 50 : 50 beamsplitter 
are related by 



giCQ-a^t+T) E 2 (t)-iE l (t-T) 

Ut) = -j, , Ht) = , (8) 

where ±cT is small difference of path length in the two arms. In Section SI and S3 of 
the SI we show that Eq. (7) may be recast as 

1 f 00 

R 34 (r r ,r s ; t) = — =■ d 2 r r d 2 r s 

c ^ ) (2ti) 2 J-oo 

[< ) (r r ,r;;o)< ) (r s ,r;,o)^' ) (r;,T0+ (9) 
w$ (r r ,T'n - t) w$ (r s ,r;,r)< } (r;,r;)] 

+ (s^r,T^-T). 

Here rj = { tj,coj} represents the set of parameters of the matter plus field incident on 
the detector; = r, s. Eq. (9) is given by the spectral and temporal overlap of the Wigner 
spectrograms of detectors W$ Eqs. (S3) (Section SI) and bare signal pathways 
R$ and (S12)-(S13) (Section S3) of SI. The detector provides a window of 
observation centered at time tj and frequency cbj,j = s, r. tj and (bj can be varied 
independently. However the temporal er^ and spectral o ] m resolutions are not 
independent and must satisfy a^l c? T > I 22 . We assume a point-size detector and 

omit all effects of spatial resolution only retaining temporal and spectral 

gating 25 ' 47 . 
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Spectral diffusion. For an electronic system coupled to a harmonic bath the matter 
correlation function obtained by the second order cumulant expansion 30 (see Section 
S2 of SI) 



F a (ti,M3,t4) = |^|V to -( tl - fe+ft -^e®-( tl ' feAA ) : 



(10) 



where (o a = co £y —co gx is the absorption frequency, ® a (fi, t 2 , t 3 , t 4 ) is the four-point 
lineshape function d> a (t lf t 2 , t 3 , t 4 ) = - t 2 ) - g a (t 3 - t 4 ) + - t 3 ) - gj,t 2 - 

h) + gdJi ~ t 4 ) ~ gMi ~ t 4 )- We shall use the overdamped Brownian oscillator 
model for the spectral density. The lineshape function then depends on two 
parameters: the reorganization energy 2 a and the fluctuation relaxation rate A a (see 
Section S2 of the SI) in the high temperature limit /c 5 r»/zA a we have 



.K 



' + A«,f-1 



(11) 



For a given magnitude of fluctuations A a , a = a, b the FWHM of the absorption 
linewidth 30 



2.355 + 1.76(A a /A a 



1 + 0.85(A a /A a ) + 0.88(A a /A«) i! 



(12) 
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